Spin-resolved impurity resonance states in electron-doped cuprate superconductors 
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With tlie aim at understanding the non-monotonic d^2_j^2-wave gap, we analyze the local elec- 
tronic structure near impurities in the electron-doped cuprate superconductors. We find that the 
local density of states near a non-magnetic impurity in the scenario of d^2_y2-wave superconduc- 
tivity with higher harmonics is qualitatively different from that obtained from the d^2_y2-wave 
superconductivity coexisting with antiferromagnetic spin density wave order. We propose that spin- 
polarized scanning tunneling microscopy measurements can distinguish the two scenarios and shed 
light on the real physical origin of a non-monotonic d^2_y2-wa.Ye gap. 
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Pairing symmetry in the electron-doped cuprate high- 
temperature superconductors has been extensively stud- 
ied experimentally and theoretically. In contrast to the 
hole-doped cuprates, where the d3.2_„2-wave pairing sym- 
metry has been generally acceptedliS,, the earlier point 
contact tunneling spectra suggested an s-wave like sym- 
metry due to the absence of zero bias conductance peak 
in the spectrum'^'"'. Recently, the phase sensitive scan- 
ning (SQUID) measurements-, nuclear magnetic reso- 
nance study^, and angle-resolved photocmission spec- 
troscopy (ARPES) experiments^'* have provided strong 
evidences that the electron-doped cuprates are also the 
(ia;2_j^2-wave superconductors. However, the functional 
form of the dj;2_y2-wave gap in electron-doped materi- 
als is a more subtle issue. A non-monotonic dj.2_j,2-wave 
gap with a maximal value in between nodal and antin- 
odal points on the Fermi surface (FS) has been measured 
in Raman experiments in NCCO^ and ARPES data on 
the leading-edge gap in Pro.89LaCeo.iiCu04-. 

Up to now, the physical origin of such non-monotonic 
da;2_j^2-wave gap is still under debate. Two kinds of the- 
oretical explanations have been put forward. One is to 
extend the superconducting (SC) gap out of the simplest 
dx2-y2-waMe via the inclusion of higher harmonics ^'-2,. 
Based on the theoretical assumption, the (ij.2_y2-wave 
pairing is caused by the attractive interaction with the 
continuum of overdamped antiferromagnetic (AF) spin 
fluctuations, which generates a maximal gap near the 
hot spots (the points along the FS separated by the 
AF move vector Qaf)- Since the hot spots in the op- 
timally doped NCCO and PCCO are located close to 
Brillouin zone diagonals, one can generally expect the 
d2;2_j,2-wave gap to be non-monotonic. The other one is 
the coexisting scenario in which the AF long-range order 
coexists with the dx2_y2-wa.ve order^iS^ii. The neutron 
scatteringi^ and transport experiments^^ have observed 
a robust AF order, which survives a broad doping region 
in the phase diagram. On the other hand, the ARPES 
measurements revealed the intriguing doping evolution 
of the FS in NCCOi^, where two inequivalent pockets 
around (7r,0) and (7r/2,7r/2) shown in the FS have been 
explained to the band folding due to the AF order^^. 
As a consequence, the resulting quasiparticle excitation 



can be gapped by both orders, and the non-monotonic 
dx2-y2-waYe gap appears naturally. 

In this paper, we argue that the local electronic struc- 
ture near impurities can provide important insight into 
the physical origin of non-monotonic dx2_y2-weLve gap. 
We calculate local density of states (LDOS) around 
a non-magnetic impurity starting from two scenarios: 
dj.2_y2--waYe superconductivity with a higher harmonic 
versus c?2;2_j,2-wave superconductivity coexisting with AF 
spin density wave (SDW) order. We find that the behav- 
ior of density of states (DOS) in both scenarios suggests 
the presence of a non-monotonic dx2_y2--wave gap. Tak- 
ing the single non-magnetic impurity into account, we 
find that in the scenario of c?;j,2„j,2-wave superconductiv- 
ity with a higher harmonic, the LDOS behaves similar to 



that shown in hole-doped cupratei 
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a single reso- 



nance state near Fermi energy appears at impurity site. 
However, due to introducing AF SDW order in the lat- 
ter scenario, the LDOS indicates a spin-resolved feature, 
i.e. two resonance states occur at impurity site with dif- 
ferent energies. For the sufficiently large SDW order, 
one spin component (spin-up or spin-down) turns out to 
be dominant, and although the DOS shows a U-shaped 
behavior, the presence of resonance states at low ener- 
gies in LDOS rules out the possibility of s-wave pairing 
symmetrjj^iiiilii^. Thus, we conclude that the different 
electronic structure near a non-magnetic impurity can 
differentiate between above scenarios and can be checked 
by the further scanning tunneling microscopy (STM) ex- 
periments. 

We start from a phenomenological superconducting 
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Hamiltonian Hsc = Z^kcrKkCk^Ckcr + Ak(cJ^|^_ki ' 
c-kj.Ck|)], where c^^ (cko-) is the fermion creation 
(destruction) operator for an electron in the state 
with wave vector k and spin projection cr, and 
^k = ^k — M with the normal state tight binding disper- 
sion Ek = —2t[cos{kx) + cos{ky)] — Ml cos{kx) cos{ky) — 
2t2[cos(2fc:i,) + cos(2fcy)] — 4t3[cos(2fca;) cos(fcjy) -I- 
cos(fc2:) cos(2/cj,)] — 4i4 cos(2fca;) cos(2fcj,) where 
{t,ti,t2,t3,U,fi) = (120,-60,34,7,20,-82) with 
the unit of meV at 0.11 dopinglS reproduce the under- 
lying FS as inferred from recent ARPES experiment*!. 
As argued above, the maximum SC gap is achieved 
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near hot spotaii^ii^, which are located much closer to 
the zone diagonal, leading to a non-monotonic behavior 
of the SC gap. A good fit of Ak to the experimental 
data is achieved via the inclusion of a higher harmonic, 
such that Ak — ^^^j^ 3 A.i[cos(ifc;j,) — cos{iky)]/2 with 
Ai ~ 5.44 meV and A3 = —2.34 meV ensures that the 
maximum of |Ak| along the FS is located at the hot 
spots. Corresponding FS and a non-monotonic gap as 
a function of the FS angle have been shown in Fig. 4 of 
Ref.20. 

By introducing a two-component Nambu spinor oper- 
ator, ^fk = (ckTi cl[^|)^, the matrix Green's function Go 
in the superconducting state can be written by 

iuJnTQ + CkT2 + AkTi 



Go(k, iw„ 



(1) 



with £'k — (^k + ^k)^^^ quasiparticle spectrum and 
Ti being the Pauli spin operator. The corresponding real- 
space Green's function is 



Go{i,j;iujn) 



N ^ 



ik-R 



■^-Go(k,ic^„), (2) 
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FIG. 1: (Color online) The dependence of LDOS on scattering 
strength f/o (a) at impurity site and (b) on the impurity's 
nearest-neighbor site. A„ and Am denote the antinodal gap 
and maximum gap at the hot spot respectively. 



where Ry = Ri — Rj with Ri being lattice vector. In 
the presence of a single-site nonmagnetic impurity of 
strength Uq located at the origin = 0, the site depen- 
dent Green's function in term of the T-matrix approach^^ 
can be obtained as 



Go{i - j;iujn) 

Go{i;iujn)T{iujri)Go{j;iUn), (3) 



where 



T{iuJn 



U0T3 



l-C/or3Go(0,0;zw„) 



(4) 



For the d-wave (with or without a higher harmonic) pair- 
ing symmetry, one can find that the local Green's func- 
tion Go{i,i;iuJn) is diagonal. As a result, the diagonal 
T-matrix reads 



T'ii,22(iw„) 



l-C/o[Go(0,0;±iw„)]i 



(5) 



where the plus (minus) sign denotes Tn (T22), giving 
rise to a particle- {ujres < 0) and hole-like {uires > 0) 
resonance state. These resonance states generate the 
sharp peaks in the LDOS only in the unitary limit 
(I uJres I /Ai < 1) where 1 = [/oi?e[Go(0, 0; ±cj^es)]ii- 

Finally, the LDOS which can be measured in the STM 
experiment is expressed as 



N{r,uj) 



with spin-resolved LDOS 



iVx(r,^) 



1 



ImGii (r, r; + iO^), 



-ImG22(r,r; 



-ui 



(6) 

(7) 
(8) 



The above equations allow a complete solution as long as 
the order-parameter relaxation can be ignored. 

In Fig.l we present the dependence of LDOS on scat- 
tering strength Uq. The LDOS with Uq = (thick 
solid line) which is equivalent to DOS in the clean sys- 
tem, shows two van Hove singularities at correspond- 
ing antinodal gap (A„) and maximum gap at the hot 
spot (Am), indicating the presence of a non-monotonic 
d2;2_j,2-wave gap, and is qualitatively consistent with 
the recent doping dependence of tunneling spectrum in 
Pr2_a;Cej;Cu04_^. In the presence of a non-magnetic 
impurity, a single resonance state appears at the im- 
purity site (Fig. 1(a)). With increasing Uq, the posi- 
tion of the resonance state shifts to positive high en- 
ergy; meanwhile the spectral weight gradually reduces 
(the LDOS in Fig. 1(a) with Uq = 50eV (thin sohd hue) 
has been amplified 500 times) and finally vanishes in the 
limit Uq — *-oc. In the LDOS on the impurity's nearest- 
neighbor site (Fig. 1(b)), there are two resonance states 
locating at the positive and negative energy with differ- 
ent spectral weight due to the particle-hole asymmetry. 
These features can be understood from the Eq. (3) and 
Eq. (5), where the correction to G{i,i;iujn) due to im- 
purity scattering reads 



5Gii(i;iLL)„ 



Uo[GQ{i;iUn)]n 
1 - Uo[Go{i;iLOn)]ii 

Uo[Go{i;iuJn)]l2 
1 - C/o[Go(i; -iLUn)] 
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(9) 



At the impurity site, the fact [Go{i;ii^n)]i2 = leads to 
a single resonance state; while on the impurity's nearest- 
neighbor site, spectral weight of both resonance states is 
nonzero, i.e. [Go(i; ia;„)]ii^i2 7^ give rise to two reso- 
nance states. Due to the same spin component at sin- 
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gle site, the spin-resolved LDOS (iV| and Ni) degener- 
ates, resulting in a degenerate single resonance state at 
the impurity site and two resonance states on the im- 
purity's nearest-neighbor site in total LDOS. These fea- 
tures, which are qualitatively similar to the (ij,2_y2-wave 
hole-doped cuprates^^, indicate that the inclusion of a 
higher harmonic in the gap function basically can not 
alter the local electronic structure near a non-magnetic 
impurity, although it generates a non-monotonic dx2_y2- 
wave gap in electron-doped cuprate superconductors. 

We now compare the above results with those re- 
sulted from the coexisting AF SDW and SC phase. 
It is convenient to introduce a 4 x 4 matrix for- 
mulation, take four-component Nambu spinor ip^ = 
(ckT,Ck-i-QT,cLk|,cLk_Qj,)^ with Q = (tt, tt) being the 
nesting vector, and then write the mean-field Hamil- 
tonian as Hsc+sdw = V'l^ (CkTsPo + MtiPq + 

Akr3Pi)^k, where rgpi = (^^^ ) , M is AF SDW 

order parameter and Ak = Ai[cos(fc:r) — cos{ky)]/2 is 
monotonic dx2-y2--wave SC gap. Note that from now on 
the wave vector k is restricted to the magnetic Brillouin 
zone (MBZ). 

The relevant matrix Green's function is obtained as 

gg^{'k,iun) = iujn ~ CkTsPo - MtiPq - AkTapi. (10) 

To solve for the resonance state in the coexisting AF 
SDW and SC phase, we define the 2x2 Green's function 
as 

Go(^,J;^c.„) = ^Y.^^'^''''' 

k 

where 

+ e''^-^'n5o]f(k,*c.„) + Mi"(k,zc.„), (12) 

with I = 1,2,3,4 denoting the left-top, right-top, 
left-bottom and right-bottom 2x2 block element of 
go(k, iuJn)- Applying T-matrix approacbi^, we can easily 
get the LDOS in the presence of a non-magnetic impurity. 

The dependence of LDOS on scattering strength Uq for 
different AF SDW order M is plotted in Fig. 2. Following 
the discussions in Ref.l9, we take the independent parti- 
cle dispersion ^k, and consider the doping dependent AF 
SDW order M . We in the following calculation choose 
the self-consistent value M — Q.lAeV at 0.11 dopingi^, 
and the decreasing M values corresponding to the doping 
increasing . 

Before considering the effect of the impurity we briefly 
review the evolution of DOS with AF SDW order in the 
SC stated. In Fig.2 the LDOS with Uq ^ Q (thick sohd 
line) is equivalent to DOS in the clean system. In the 
limit M = OeV, as seen in hole-doped cuprates^, the 
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FIG. 2: (Color online) The dependence of LDOS on scatter- 
ing strength Uo for different AF SDW order M. LDOS at 
impurity site: (a)Af = OeV, (b)M = 0.05eV, and (c)M = 
0.14eV; and LDOS on the impurity's nearest-neighbor site: 
{d)M = OeV, (e)M = 0.05eV, and (f)M = 0.146^. 



DOS (thick solid line in Fig. 2(a)) at low energies be- 
haves to be V-shaped like with a monotonic d^2_j^2-wave 
SC gap, and a coherent peak locates at the maximal gap 
edge. After introducing AF SDW order, another coher- 
ent peak appears at the energy less than the maximum 
gap (thick solid line in Fig. 2(b)). Thus, a non-monotonic 
drc2-y2-wave SC gap occurs in coexisting AF SDW and 
SC state. In particular, the DOS shows the U-shaped 
behavior at sufficiently large SDW order M = 0.14ey 
with doping x = 0.11^^- (thick solid line in Fig. 2(c)), 
which has been observed in earlier point contact tunnel- 
ing spectraM. These unusual evolutions of DOS with AF 
SDW order are qualitatively similar to the doping evolu- 
tion of DOS'^^, and have been explained as the result of 
the coexisting AF SDW and SC state"'20. 

We proceed to analyze the dependence of LDOS on 
scattering strength Uq near a non-magnetic impurity. 
For the limit AF SDW order M = OeV, we show 
that the LDOS at the impurity site (Fig. 2(a)) and on 
the impurity's nearest-neighbor site (Fig. 2(d)) are sim- 
ilar to that obtained from the dj.2_j^2-wave hole-doped 
superconductors^®. With increasing AF SDW order, the 
LDOS at the impurity site are qualitatively different. In 
Fig. 2(b), it is clearly shown that two resonance states at 
the impurity site occur near the Fermi energy, which in- 
dicates that the degenerate single resonance state with 
M — has separated due to the presence of AF SDW or- 
der. In principle, the spin-resolved LDOS should give rise 
to multiple resonance states on the impurity's nearest- 
neighbor site, though they are actually not easy to be re- 
solved in Fig. 2(e) because of the resonance states crossing 
each other near Fermi energy. At sufficiently large SDW 
order M = 0.14ey, one resonance state at the impurity 
site (Fig. 2(c)) exists and shifts towards the gap edge, the 
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FIG. 3: (Color online) The spin-resolved LDOS A^^ (dashed 
line) and (solid line) at impurity site without scattering 
Uo = OeV for AF SDW gap (a)M = OeV, (b)M = O-OSeF ,and 
(c)M = O.UeV; and with Uo = 0.5eV for (d)A/ = OeV, 
(e)M = 0.05c V'.and (f)M = O.UeV. 



other one is barely visible due to the vanishing spectral 
weight. For the better understanding of such important 
features, the spin-resolved LDOS {N-^ and Ni) at impu- 
rity site without scattering Uq = OeV and with scattering 
strength Uq = 0.5eV are shown in Fig. 3 for increasing 
M. When M = 0, the spin- resolved iV| and Ni with 
Uo — OeV degenerate (Fig. 3a), the resulting resonance 
state in spin-resolved LDOS Ni and Ni with Uq = 0.5eV 
are located at the same resonance energy, leading to a 
degenerate single resonance state (Fig. 3d). With M in- 
creasing, the LDOS is not equal to Ni (Fig. 3b), thus 
the degenerate spin-resolved LDOS separates, leading to 
two resonance states at impurity site with different spec- 
tral weight (Fig.3e). At sufficiently large SDW order M, 



the LDOS Ni is dominant over (Fig. 3c), thus one sin- 
gle resonance state from spin-up component exists and 
shifts towards the gap edge due to the U-shaped DOS 
(Fig.3f), and the other one from spin-down component is 
barely visible due to the vanishing spectral weight (solid 
line in Fig. 3c). Therefore due to the different spin com- 
ponents induced by the presence of AF SDW gap, an 
existing impurity will be polarized by a local net spin-up 
or spin-down, which leads to the splitting of the LDOS. 
These unique features do not appear in the scenario of 
rf2.2_j^2-wave superconductivity with a higher harmonic 
and should be detected by the spin-polarized STM mea- 
surement. 

In summary, we analyze the LDOS around a non- 
magnetic impurity in electron-doped cuprate supercon- 
ductors starting from two different scenarios: d^^-y^- 
wave superconductivity with a higher harmonic versus 
(i2.2_j^2-wave superconductivity coexisting with AF SDW 
order. We find that in both cases the DOS indicates 
the presence of non-monotonic dj.2„j,2-wave gap, quali- 
tatively consistent with the recent tunneling spectrum 
measurement in Pr2_a;Cea;Cu04_5^^, therefore both of 
them have been thought to be the possible physical ori- 
gins of the non-monotonic c?2,2_j,2-wave gap. We also find 
that the inclusion of a higher harmonic basically doesn't 
alter the local electronic structure near a non-magnetic 
impurity; in contrast, with introducing AF SDW order, 
the LDOS presents spin-resolved feature, i.e. a degener- 
ate single resonance state at the impurity site separates 
into two resonance states due to the different spin com- 
ponent induced by the presence of AF SDW order. Thus 
we strongly suggest that the future spin-polarized STM 
measurements should be performed to differentiate two 
above scenarios and shed light on the real physical ori- 
gin of non-monotonic dj.2_y2-wavc gap in electron-doped 
cuprate superconductors. 
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